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摘要 : 本 文通 过 对 肿瘤 与 生物 组 织 的 物理 模型 简化 ， 应 用 生物 传 热 方程 ， 建 立 肿 瘤 热 疗 的 数学 模型 。 运 用 数值 模拟 预演 了 热 疗 
过 程 ， 通 过 控制 探 针 的 开关 ， 控 制 了 肿瘤 中 心 温度 在 54C 以 下 ， 肿 瘤 边 缘 温 度 在 43~45C 之 间 变 化 ， 分 析 了 肿瘤 内 部 的 温度 分 
布 ， 提 出 了 基于 温度 控制 的 热 疗 方案 ， 为 临床 肿瘤 热 疗 的 温度 控制 提供 理论 依据 。 
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A Numerical Study on Temperature Research for Microwave Thermal 


Therapy of Tumor 
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Abstract: This article simplifies the physical model of tumor and biological tissue, applies the bio heat transfer 
equation to build a mathematic model of tumor hyperthermia. A numerical method is adopted to previews the tumor 
treatment process, the tumor center temperature under 54'C, the tumor edge changes between 43~45'C by 
controlling the switch of the antenna, the temperature distribution within the tumor is also analyzed, put forward the 
heat treatment scheme based on temperature control. The article provides theory foundations for temperature control 
of clinical tumor thermotherapy. 
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0 引 言 ' 


热 疗 法 是 近年 来 迅速 兴起 的 一 种 有 效 治疗 肿瘤 。” 度 场 预 示 。 程 志 刚 等 中 根据 体 模 实验 所 得 到 的 SAR 
的 手段 。 热 疗法 具有 无 创 或 微 创 的 特点 ， 无 毒 无 副 。 公式， 利用 有 限 元 软件 ANSYS 模拟 出 915 MHz 相 
作用 ， 因 而 治疗 过 程 中 病人 的 痛苦 较 小 。 热 疗法 不 ”入 式 探 针 的 微波 热 场 ， 同 时 考虑 到 微波 加 热 体 模 后 
但 能 直接 杀 死 癌 细胞 ， 而 且 可 作为 放疗 、 化 疗 的 辅 。 “变形 造成 的 部 分 参数 的 变化 ， 对 模拟 中 的 相关 参数 
助 治疗 方法 。 进行 了 修改 ， 从 而 得 到 了 与 实际 测量 较为 接近 的 温 

长 期 实践 证 明 ， 微 波 热 疗 优 于 其 他 各 种 物理 加 升 曲线 。T P Ryan 采用 三 维 有 限 元 模型 预测 了 腰痛 
热 方法 。 介 入 微波 热 凝固 疗法 是 欧美 等 国 进行 热 疗 ”治疗 中 少 椎 组 织 附近 电极 产生 的 电场 和 热 沉 积 作用 
的 主要 方法 , 对 肿瘤 有 较 好 的 疗效 。 近 年 来 已 开发 出 “的 能 量 分 布 ， 发 现 模拟 结果 与 活体 动物 以 及 人 体 标 
应 用 微波 加 热 深部 肿瘤 的 设备 中 , 实现 了 对 体内 比 吸 。 ”本 实验 的 结果 一 致 。 

收 率 SAR (specific absorption rate) 分 布 的 控制 而 达到 在 上 述 作者 的 研究 基础 上 ， 本 文 以 肝脏 肿瘤 为 
对 不 同 深度 的 肿瘤 进行 加 热 。 朱 光明 等 中 以 前 列 腺 癌 ”对 象 ， 采 用 有 限 元 方法 ， 分 析 热 疗 过 程 的 温度 场 ， 
热 疗 为 研究 考察 对 象 ， 讨 论 介入 式微 波 热 疗 对 微波 “从 而 提出 了 基于 温度 控制 的 热 疗 方案 。 

热 疗 过 程 中 的 温度 场 、 热 损伤 分 布 以 及 血液 灌注 率 a 

的 变化 进行 了 数值 模拟 ， 得 到 了 较为 准确 的 热 疗 温 。 | 微波 八 疗 贷 型 
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1.1 微波 热 疗 物理 模型 

以 肝 肿 瘤 热 疗 为 研究 对 象 。 在 临床 研究 中 ， 对 
于 D<2 cm 的 肝 肿 瘤 ， 采 用 单一 热源 ， 一 次 凝固 即 可 
完成 治疗 。 本 文 假定 肿瘤 为 直径 D=1.5 cm 的 球体 ， 
采用 单一 探 针 加 热 。 

微波 热 凝 固 肿瘤 热 疗 系统 模型 如 图 1 所 示 ( 训 
面 图 ), 取 915 MHz 圆柱 形 微波 探 针 ,半径 Ri=0.001 
m， 前 端 微波 有 效 发 射 长 度 H=0.008 m; 计算 区 域 取 
以 肿瘤 中 心 为 球 心 、 半 径 R=0.045 m 的 球体 ， 肿 瘤 
周围 计算 区 域 为 正常 肝 组 织 ， 探 针 有 效 发 射 端 插入 


肿瘤 中 心 ， 圆 柱 体 中 心 与 肿瘤 中 心 重合 。 
日 
肿 交 ® 
所 月 如 
R R 
探 和 有 这 区 


图 1 ”肝脏 肿 瘤 微 波 热 疗 示 意图 


Figure.1 Schematic diagram of the liver tumor microwave 


hyperthermia 


1.2 控制 方程 及 定 解 条 件 
本 文采 用 生物 传 热 中 最 为 广泛 应 ) 
星 ， 如 式 (1) 所 示 。 


式 中 ，p 为 组 织 或 肿瘤 的 密度 ，(kg/m’); c 为 组 织 或 
肿瘤 的 比热容 ，(J/kg*K); 7 为 计算 时 间 ，(S); 7 为 温 
度 ，(K); 4 为 导热 系数 ，(W/m*K) ; 1 为 血液 治 注 
率 ，(kg/m*s); CG; 为 血液 比热容 ，(J/kg 民 ); 为 血 
液 温 度 ，('C)，g 为 组 织 的 代谢 产 热 率 ， (W/m’); gq; 
为 单位 组 织 所 吸收 的 微波 辐射 能 ，(W/m)。 

计算 区 域 边界 温度 取 为 定 值 ， 热 疗 前 的 温度 场 
为 初始 温度 场 。 由 于 加 热 过 程 产生 的 热 场 是 绕 探 针 
轴 对 称 的 ， 本 问题 可 以 采用 二 维 数值 模拟 。 
1.3 物性 参数 

血液 灌注 率 瑟 定义 为 单位 时 间 单 位 体积 组 织 内 
的 血液 流量 ， 是 活体 组 织 区 别 于 非 生物 材料 的 重要 
特性 。 血 液 灌 注 率 按 以 下 所 示 式 (2)、 式 (3) 计 算 9。 


的 Pennes 方 
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(0-45.0 ) 

0.45+3.55*e 120 7T<45.0 (2) 
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0 T>60 
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4.8 

0.833- 亿 370 370<7<42.0 G3) 
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比 吸收 率 SAR 是 用 来 衡量 电磁 辐射 对 人 体 的 影 
响 的 重要 指标 .SAR 定义 在 单位 时 间 内 每 千克 的 物质 
所 吸收 的 电磁 能 量 。 计 算 中 比 吸收 率 SAR 采 用 915 
MHz 圆 柱 形 微 波 探 针 实 验 数 据 门 。 
由 于 肿瘤 组 织 代 谢 旺 盛 ， 因 此 肿瘤 与 正常 组 
代谢 产 热 率 差 距 较 大 ， 肿 瘤 与 肝 组 织 的 代谢 产 热 
别 取 gwwwwo=29000 W/o 和 gyjwo=450 W/m 名。 

其 他 物性 参数 取 常 数 ， 具 体 取 值 如 表 1 所 示 纪 。 


表 1 物性 参数 
Table 1 The physical parameters 


让 AS 
孙 站 


p 4 c qm 

kglm W/mK) kgeK) WA 
组 织 1020 0.465 3840 450 
肿瘤 1000 0.642 3500 29000 
血液 1060 一 一 3770 


2 数值 模拟 


2.1 热 疗 前 的 初始 温度 场 

把 控制 方程 (D) 式 中 等 式 左 侧 的 非 稳 态 项 和 等 式 
右 侧 第 四 项 分 别 置 为 零 ， 就 可 用 于 计算 热 疗 前 的 温 
度 场 ， 即 模拟 热 疗 过 程 的 初始 温度 场 。 模 拟 所 用 的 
网 格 划分 及 模拟 得 到 的 初始 温度 场 如 图 2 所 示 ， 肿 瘤 
中 心 最 高 温度 为 到 41.29C ， 这 是 由 于 肿瘤 相对 于 正 
常 组 织 而 言 代 谢 更 旺盛 所 致 。 


Temperature initial field 


图 2 模型 网 格 划 分 与 初始 温度 场 
Fig.2 The mesh and the initial temperature field of the model 


2.2 探 针 加 热 后 的 温度 场 


在 初始 温度 场 分 布 的 基础 上 ， 加 入 微波 热源 。 
持续 加 热 100 s 后 的 温度 场 如 图 3 中 左 图 所 示 , 可 以 看 
出 ， 此 种 探 针 形成 了 沿 探 针 前 后 方向 不 对 称 的 椭 球 
形 温度 场 ， 相 同 半径 位 置 处 垂直 于 探 针 方向 的 温度 
要 明显 低 于 沿 探 针 方向 的 温度 。 为 此 本 文选 取 肿 瘤 
中 心 点 1 作为 肿瘤 中 心 温 度 监测 点 ， 选 取 肿 瘤 边 缘 外 
1 mm 的 点 2、 点 3、 点 4 作为 肿瘤 灭 活 监测 点 ， 具 体 布 
置 如 图 3 中 右 图 所 示 。 


图 3 100 s 后 的 温度 场 和 监测 点 、 路 径 示 意图 
Fig.3 The temperature field for heating 100 s by antenna and 


schematic diagram of measuring point and path 


图 4 所 示 为 微波 持续 加 热 过 程 中 监测 点 温度 随 
时 间 变 化 曲线 ， 图 5 所 示 为 加 热 100 s 后 沿 半径 方向 
(图 3 中 右 图 所 示 路 径 位 置 ) 的 温度 分 布 曲线 。 


40 60 
时 间 time/s 
图 4 监测 点 的 温度 随时 间 变 化 曲线 


Fig.4 The time histories of the monitoring points’ temperature 


温度 temp/'C 


i | i ss 
0.005 0.01 0.015 0.02 0.025 0.03 0.035 0. 
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图 5 沿 半径 的 温度 分 布 曲线 


Fig.3 The temperature distribution curve along radius 


从 图 4 温度 曲线 可 以 看 出 ,， 当 探 针 持续 加 热 100s 
后 ， 肿 瘤 中 心 点 1 温度 已 达 137C ， 此 温度 会 使 肿瘤 
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心 失 水 炭化 ， 影 响 之 后 的 热 疗 进程 ， 同 时 三 个 肿 
瘤 灭 活 监 测 点 的 最 低温 度 值 也 已 达 92'C， 严 重伤 害 
到 周围 组 织 。 从 图 5 温度 分 布 曲线 中 可 以 看 出 沿 着 半 
径 方向 温度 逐渐 减低 ， 在 半径 为 0.025 m 时 温度 趋 于 
定 值 37.S'C 。 


3 基于 温度 控制 的 热 疗 模 拟 


在 微波 热 疗 过 程 中 ， 温 度 控 制 是 通过 实时 监测 
来 实现 的 ， 然 而 在 现 有 条 件 下 ， 有 损 测量 会 给 病人 
增加 痛苦 ; 无损 测量 的 精度 暂时 无 法 满足 治疗 过 程 
的 要 求 。 针 对 以 上 缺陷 ， 本 文 运 用 数值 模拟 预演 
热 疗 过 程 ， 提 出 基于 温度 控制 的 热 疗 方案 ， 实 现 既 
能 杀 死 整个 肿瘤 组 织 ， 又 能 保证 周围 正常 组 织 不 受 
过 度 损伤 。 
3.1 温度 控制 原则 

为 了 获得 一 个 理想 的 肿瘤 热 疗效 果 ， 加 热 时 必 
须 满 足以 下 条 件 : 在 限定 的 时 间 内 将 肿瘤 组 织 加 
热 到 43'C 以 上 ， 并 维持 一 定时 间 ; @ 在 肿瘤 组 织 内 
不 能 存在 41'C 以 下 的 冷 点 ，@ 肿 瘤 周围 组 织 温度 必 
须 保持 在 45C 以 下 外。 
3.2 控制 方案 

从 图 4 可 以 看 出 当 探 针 加 热 14 s 时 ， 点 4 的 温度 为 
45C， 点 2、 点 3 的 温度 为 S1C 左 右 ， 点 1 的 温度 为 
54C， 对 应 的 温度 场 分 布 如 图 6 中 左 图 所 示 。 因 此 ， 
为 了 保证 热 疗 的 效果 ， 满 足 温度 控制 原则 ， 必 须 控 
制 使 点 4 的 温度 满足 热 疗 的 要 求 ， 此 时 点 2、 点 3 的 温 
度 也 不 至 于 太 高 ， 不 会 过 度 伤害 正常 组 织 。 


i 
py ES 37.5 
\ 38.4209 
. 39.3418 
40.2627 
41.1836 
42.1046 
43.0255 
43.9464 
\ / 44.8673 
醒 了 45.7882 
ure field Temperature field 


图 6 热 疗 方案 进行 14 s、200 s 时 的 温度 场 
Fig.6 Temperature field after 14 s and 200 s in heat treatment 


plan 


按照 温度 控制 原则 ， 加 热 14 s 后 关闭 微波 探 针 
停止 加 热 ， 让 热量 从 中 心 高 温 处 逐渐 向 外 围 低温 处 
传导 ， 此 时 由 于 血液 的 影响 ， 将 带 走 中 心 热量 ， 从 
而 使 肿瘤 与 组 织 的 温度 降低 。 到 200 s 时 ， 对 应 的 温 
度 场 如 图 6 中 右 图 所 示 ， 点 2 的 温度 下 降 到 43.8°C， 
点 3 温度 下 降 到 44.1C， 点 4 下 降 到 43C， 点 1 下 降 到 
45.8'C， 满 足 温 度 控制 的 要 求 ， 能 够 杀 死 肿瘤 而 不 


损坏 正常 组 织 。 
3.3 热 疗 模拟 


制定 了 基于 温度 控 


通过 对 模拟 过 程 中 相应 的 监测 点 的 观测 ， 本 文 


二 
上 


dt 


上 的 热 疗 模拟 ， 既 能 控制 肿瘤 中 


心 点 温度 不 会 过 高 


而 炭化 ， 也 能 控制 肿瘤 边缘 温度 ， 


达到 彻底 消灭 肿瘤 


不 损伤 正常 组 织 的 要 求 。 对 温 


度 控制 ， 是 通过 控 


探 针 ， 以 此 控制 整 
程 温 度 变 化 曲线 如 
波 探 针 的 加 热 时 间 


波 探 针 ， 肿 瘤 边缘 监测 点 4 从 初始 温度 37.5'C 升 高 至 


45C 后 关闭 探 针 ,， 待 监测 点 4 温度 降 至 43C 后 再 打开 


制 探 针 的 开关 来 实现 的 。 接 通 微 


个 热 疗 过 程 的 温度 。 热 疗 模拟 过 
图 7 所 示 ， 可 以 看 出 ， 首 次 接 通 微 
为 14 s， 关 闭 探 针 的 冷却 时 间 为 


186 s， 后 续 加 热 周 期 约 为 4 s 左 右 , 冷却 周期 约 为 152 


s 左 右 。 
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7 监测 点 基于 温度 控制 随时 间 变 化 曲线 


Fig.7 The time histories based on the monitoring points 


4 结论 


本 文 提出 了 基 


于 温度 控制 的 热 疗 方案 ， 通 过 数 


值 模拟 预演 了 915 MHz 的 单 点 介入 式微 波 热源 杀 死 
肿瘤 的 热 疗 过 程 ， 数 值 模拟 结果 表明 ， 可 设 定 初期 
微波 探 针 加 热 时 间 14 s， 之 后 冷却 至 200 s， 随 后 呈 周 
期 性 加 热 与 冷却 ， 加 热 时 间 约 为 4 s， 冷 却 时 间 约 为 
152 s， 可 控制 肿瘤 中 心 最 高 温度 在 54'C 以 下 ， 肿 瘤 
边缘 温度 在 43~45C 之 间 变 化 ， 达 到 控制 肿瘤 中 心 温 
度 不 致 过 高 、 同 时 杀 死 肿瘤 又 不 会 过 度 损伤 周围 组 


织 的 目的 。 
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